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FOREWORD 


This  report  represents  one  phase  of  a  general  study  of  the 
Interaction  of  Ions  with  fluid  flows  under  Contract  DA-31-124- 
ARO-D-246,  U.  S.  Army  Research  Office  -  Durham,  with  Dr.  Henry 
R.  Velkoff  serving  as  principal  Investigator. 

The  authors  wish  to  thank  T.  Chuang,  X.  Dlller,  D.  Smith, 
I.  Anderson,  and  M.  Wlnlass  for  their  assistance  In  various 
phases  of  the  work  reported  herein. 
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ABSTRACT 


Prletlon  factors,  velocity  distributions,  pressure  gra¬ 
dients,  and  density  measurements  were  made  on  a  5/8-lnch  x 
5-inch  rectangular  flow  channel  with  air  flow.  A  plane  of 
corona  wires  along  the  midplane  of  the  channel  was  the  source 
of  Ions  which  drifted  across  the  air  flow  to  the  electrically 
grounded  walls  of  the  channel.  Experimental  data  were  collected 
to  determine  the  Influence  of  the  Ion  flow  on  the  air  flow  and 
the  extent  of  the  lon-fluld  coupling. 

The  electrodynamic  and  hydrodynamic  equations  were  set  up 
and  several  approximate  solutions  were  obtained  based  on  the 
postulation  of  a  flow  induced  body  force  and  the  assumption  of 
fully  established,  one-dlmenslonal  flow.  The  magnitude  of  the 
Ion  effects  was  found  to  be  significantly  dependent  on  the  Ion 
density  distribution  within  the  flow  channel. 


ill 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I  INTRODUCTION  1 

II  APPROXIMATE  THEORY  2 

Hydrodynamic  Approach  2 

Kinetic  Theory  Approach  12 

III  EXPERIMENTAL  l£ 

Description  of  Flow  Channels  l£ 

Electrical  Characteristics  of  the  Test 
Apparatus  17 

Effect  of  Ions  on  Friction  Factor  28 

Effect  of  Ions  on  Transverse  Pressure  Gradient  37 

Effect  of  Ions  on  Velocity  Distribution  38 

Effect  of  Flow  on  Current  Density  43 

Flow  Visualization  Attempts  51 

Supplementary  Tests  55 

IV  CONCLUSIONS  AND  RECOMMENDATIONS  58 

REFERENCES  60 

APPENDICES 

I  HYDRODYNAMIC  AND  ELECTRODYNAMIC  EQUATIONS  63 

II  APPROXIMATE  ANALYSIS  SATISFYING  THE  IRROTATIONAL 

ELECTRIC  FIELD  CONDITION  65 

III  EFFECT.  OF  NEARBY  SURFACES  ON  CORONA  67 


iv 


LIST  OF  FIGURES 


PAGE 


Friction  factor  'times  Reynolds  number  versus  charge 

number  for  two  theoretical  models  9 

Nondimensional  velocity  profile  showing  the  effect 

of  charge  number,  based  on  a  theoretical  model  13 

Schematic  drawing  of  the  flow  channel  18 

Cross  sectional  drawing  of  the  flow  channel  19 

Electrical  circuit  of  the  flow  channel  test  apparatus  20 

Photograph  of  the  flow  channel  and  instrumentation  21 

Photograph  of  the  flow  channel  22 

Corona  current  versus  volts  applied  to  wire  showing 

the  variation  in  characteristic  among  the  wires  23 

Corona  current  versus  volts  applied,  for  wire  Ko.  5, 

showing  the  different  characteristic  for  increasing 

and  decreasing  voltages  25 

Corona  current  versus  volts  applied,  for  wire  No.  5, 
showing  that  Reynolds  number  has  no  effect  26 

Corona  current  versus  voltage  applied,  for  wire  No,  5, 
showing  the  effect  of  a  period  of  time  on  the  corona 
characteristic  27 

Nondimensional  static  pressure  versus  distance  along 

duct  for  a  Reynolds  number  cf  approximately  140,  with 

charge  number  as  parameter  29 

Nondimensional  static  pressure  versus  distance  along 

duct  for  a  Reynolds  number  of  approximately  265,  with 

charge  number  as  parameter  30 

Nondimensional  static  pressure  versus  distance  along 

duct  for  Reynolds  number  of  approximately  400,  with 

charge  number  as  parameter  31 

Nondimensional  static  pressure  versus  distance  along 

duct  for  a  Reynolds  number  of  approximately  540,  with 

charge  number  as  parameter  32 


v 


LIST  OF  FIGURES  (Continued) 


] 

-j 
3 1 

0 

a 

n 


FIG.  NO.  PAGE 


16 

Nondimensional  static  pressure  versus  distance  along 
duct  for  a  Reynolds  number  of  approximately  1100,  with 
charge  number  as  parameter 

33 

17 

Log  friction  factor  versus  log  Reynolds  number,  with 
charge  number  as  parameter 

34 

18 

Friction  factor  times  Reynolds  number  versus  total 
corona  current  shoving  comparison  between  experiment 
and  theoretical  models 

35 

19 

Friction  factor  times  Reynolds  number  versus  charge 
number,  showing  comparison  between  experiment  and  a 
theoretical  model 

36 

20 

Electric  pressure  versus  nondimens ional  transverse 
distance,  tj,  for  location  A  and  C 

39 

21 

Electric  pressure  versus  nondimens  ional  transverse 
distance,  tj,  for  location  A  for  two  different  probe 
orientations 

40 

22 

Electric  pressure  versus  nondimens ional  transverse 
distance,  tj ,  for  location  B 

4l 

23 

Difference  between  pressure  at  wall  and  probe  pressure 
at  t)  *  ,5  versus  corona  current  density 

42 

24 

Nondimens  ional  velocity  profile,  comparing  experimental 
data  with  a  theoretical  model 

44 

25 

Nondimens ional  velocity  profiles  for  Reynolds  number 
of  265,  for  various  charge  numbers 

45 

26 

Nondimens ional.  velocity  profiles  for  Reynolds  number 
of  540,  for  various  charge  numbers 

46 

27 

Nondimensional  velocity  profiles  for  Reynolds  number 
of  1095 »  for  various  charge  numbers 

47 

28 

Nondimensional  velocity  profile  at  location  B 

48 

29 

Nondimensional  velocity  profile  at  location  C 

49 

30 

Fhotograph  of  the  flow  channel  with  one  wall  removed, 
showing  the  segmented  electrod-. 

52 

D 


LIST  OF  FIGURES  (Continued) 


FIG.  NO.  PAGE 


31 

Current  density  measured  at  lower  wall  versus  axial 
distance  with  Reynolds  number  as  parameter 

53 

32 

Current  density  at  lower  wall  versus  Reynolds  number, 
for  two  different  values  of  axial  distance  Z 

54 

33 

Static  pressure  on  upper  wall  versus  axial  distance 

Z/L  showing  effect  of  discontinuous  electric  field 
along  the  duct 

56 

34 

Static  pressure  on  upper  wall  versus  axial  distance  Z 
showing  the  effect  of  number  of  corona  wires 

57 

35 

Corona  current  versus  voltage  applied,  showing  the 
effect  of  a  nearby  surface 

68 

36 

Corona  current  versus  voltage  applied  to  a  needle 
point,  showing  the  effect  of  geometry 

70 

37 

Corona  current  versus  voltage  applied  to  a  needle 
point,  showing  the  effect  of  a  superimposed  air 
stream 

71 

vii 


LIST  OF  SYMBOLS 


D 


dio 

Ec 

% 

Ez 

-♦ 

F 

H 

H 

Jy 

Jz 

— » 

J 

k 

K 

L 

m* 


NPC 

Ni 

Na 

N3 

nRE 

n0 

P 

P* 

R 


Diffusion  coefficient 

Collision  diameter 

Critical,  field  strength 

Field  strength  in  y  direction 

Field  strength  in  z  direction 

Body  force 

Half  height  of  duct 

Magnetic  field  strength 

Current  density  in  y  direction 

Current  density  in  z  direction 

Current  density 

Boltzmann's  constant 

Mobility 

Length  of  duct 

Reduced  mass 

Charge  number 

Charge  number 

Charge  number 

Charge  number 

Reynold's  number 

Density  of  neutral  molecules 

Pressure 

Dimensionless  pressure,  P/fow^ 
Gas  constant 


viii 


LIST  OF  SYMBOLS  (Continued) 


Temperature 

Velocity  component  in  y  direction 

Velocity 

Volts 

Velocity  component  in  z  direction 
Mean  velocity 

Dimensionless  velocity,  w/wB 
Distance  coordinate  across  the  flow 
Distance  coordinate  along  the  flow 
Z/L 

Permittivity 

Z/L 

Viscosity 

Dielectric  permeability 
Y/H 

Fluid  density 
Ion  density 
Mean  ion  density 
Electrical  conductivity 


SECTION  I 


INTRODUCTION 

The  investigation  presented  in  this  report  is  concerned 
with  the  flow  of  slightly  ionized  air  in  a  rectangualr  channel 
in  the  presence  of  an  applied  electric  field.  The  work  pre¬ 
sented  is  an  extension  of  work  in  which  flow  of  ionized  air 
in  a  round  tube  was  studied  and  reported  . 1,2,3  The  prior  work 
utilized  gas  flowing  along  a  1-1/4-inch  diameter  tube  in  which 
a  fine  wire  was  located  concentrically.  A  voltage  was  applied 
between  the  wire  and  the  wall  of  sufficient  intensity  to  pro¬ 
duce  a  corona  discharge  around  the  wire  and  consequently  pro¬ 
duce  unipolar  ions  in  the  gas.  The  results  of  this  prior  work 
indicated  large  changes  in  pressure  drop,  velocity  profile,  and 
heat  transfer  occurred  when  slgnlflcnnt  ion-currents  flowed 
transverse  to  the  flow.  The  analysis  reported  in  Reference  3 
postulated  the  existence  of  an  axial  body  force  as  a  limiting 
case  and  the  results  of  the  analysis  showed  order  of  magnitude 
agreement  with  the  test  data.  From  the  analytical  work  it 
could  be  concluded  that  the  ion-flow  coupling  mechanism  could 
possibly  be  active  in  parallel  plate  channel  flow  and  in  ex¬ 
ternal  boundary  layers  as  well.  Consequently,  studies  of  the 
ion-flow  coupling  effects  on  boundary  layers  was  extended  into 
parallel-plate  channel  flow  and  external  flow.  This  report 
covers  the  studies  conducted  on  the  parallel  plate  case,  which 
were  aimed  at  determining  whether  the  trends  predicted  by  the 
theory  for  the  flat  channel  would  be  found  to  establish  the 
validity  of  the  theory,  and  above  all  to  secure  a  better 
understanding  of  the  basic  mechanisms  underlying  the  ion-flow 
coupling  phenomenon. 

In  the  present  work,  theoretical  analyses  and  experimental 
tests  were  carried  on  using  air  at  atmospheric  pressure  and 
temperature  with  a  corona  discharge  as  the  source  of  ionization. 
Experimental  results  show  that  the  friction  factor  of  the  fluid 
flowing  in  the  duct  is  Increased  by  the  ion  coupling  by  a  fac¬ 
tor  of  approximately  two,  within  the  range  of  tests  performed. 
Approximate  solutions  have  been  developed  which  show  similar 
trends  as  the  experimental  results. 


SECTION  II 
APPROXIMATE  THEORY 

Hydrodynamic  Add  rote  h 

The  oat*  of  th*  flow  b*tw**n  two  lnflnlt*  parall*l  flat 
pianos  la  considered.  This  nodtl  approximates  flow  In  a  rec¬ 
tangular  channel  with  a  large  aspect  ratio.  Ions  of  one  sign 
are  assumed  to  be  emitted  from  a  midplane  and  travel  across 
the  flow  area  to  the  flat  plates. 


flat  plate  electrodes 
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ION  SOURCE  AT  MIDPLANE 


In  the  hydrodynamic  approach,  the  Ions  existing  In  the 
fluid  give  rise  to  a  body  force  per  unit  volume,  pcE,  acting 
on  the  fluid.  The  following  assumptions  are  also  used: 

1.  The  fluid  is  incompressible. 

2.  Gravity  forces  and  magnetic  forces  are  negligible. 

3.  Ionization  is  produced  only  on  the  midplane,  and 
therefore,  only  ions  of  one  sign  exist  within  the 
bulk  of  the  fluid. 

4.  Charge  diffusion  is  negligible. 

5.  Steady-state  conditions  exist  at  all  points  in  the 
flow  field. 

6.  Non-uniform  electric  field  terms  are  negligible. 

7.  Conductivity  is  small  and,  therefore,  all  current  is 
caused  by  ion  motion. 

8.  The  ions  have  a  constant  mobility. 
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The  complete  set  of  hydrodynamic  and  electrodynamic 
equations  for  Cartesian  coordinates  listed  In  Appendix  I 
then  reduce  to 


p  (u  in  ♦  w 

\  3y 


-  i£  + 

3Z 
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/  a2w  +  a2w\ 
v  ay2  Jzz  / 


♦  PCEZ 
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*y  'ay2  j*2  * 


pcEz 


3u  +  3w 
ay  az 


0  , 


PgU  ^  PgKEy  ■  Jy  | 

PC"  *  pcEEz  m  * 

Mi  ♦  Lk  .  o  , 

ay  az 

ilz  -  IE*  -  0  , 

az  3y 

3£Z  +  igy  .  £c 

3z  3y  c  * 


(1) 

(2) 

(3) 

<*> 

(5) 

(6) 

(7) 

(8) 


Up  to  this  point  In  the  analysis,  gradients  In  the  flow 
direction  and  fluid  motion  In  the  transverse  direction  are 
Included.  The  assumption  is  now  made  that  the  velocity  and 
electric  fields  are  "fully  established",  i.e.,  all  gradients  In 
the  z  direction  except  ap/3z  are  zero  and  u  ■  0.  Then  Eq.  (6) 
merely  states  that  Jy  Is  a  constant,  and  Jz  appears  In  Eq.  (5) 
only.  The  solution  of  Eq.  (1)  for  the  velocity  field  then  re¬ 
quires  a  specification  of  Jz. 

A  limiting  solution  may  be  obtained  by  taking  Jz  ■  0.  It 
Is  seen  then  that  by  combining  Eq.  (5)  into  Eo.  (7)  that 

l  -  0  ,  (9) 

K  3y 

which  is  satisfied  only  for  w  ■  0  under  the  assumptions  made. 
(See  Appendix  II  where  this  difficulty  is  removed,  but  then 
Eq.  (8)  is  not  satisfied.)  The  following  solution  then  will  be 
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a  limiting  solution  for  ths  flow  vtloclty  approaching  zero,  or 
can  be  considered  as  an  approximate  solution  for  finite  velocity 
which  does  not  satisfy  the  lrrotatlonal  field  condition  of  Eq. (7). 
Eqs.  (1)  through  (8)  become: 


If  -  -0  . 

(10) 

If  -  »cEy  -  o  . 

ill) 

Pc^Ey  ■  Jy  » 

(12) 

w  ♦  KEr  -  0  , 

(13) 

ajy  - 

-  0  , 

*y 

(14) 

JEg  _ 

r-5-  •  0  . 

Jy 

(15) 

a  Ey  Pc 

>y  "  r  • 

(16) 

Prom  Eq.  (14),  Jy  is  a  constant.  Then  substituting  pc  from 
Eq. (12)  Into  Eq.  (16),  it  Is  seen  that  pc  and  Ey  are  uncoupled 
from  the  velocity  field. 


Substituting  oc  from  Eq.  (12)  Into  Eq.  (16)  yields 

1  JEy_  . 

2  >y  Ke 

whose  solution  Is 


(17) 


Ey2  -  ♦  C1  (18) 

Kc 

where  Cx  Is  a  constant.  This  relation  has  also  been  developed 
elsewhere  (4). 

The  boundary  condition  on  Ey  specifies  the  value  of  the 
constant  Cx*  If  the  Ionization  is  described  as  "space  charge 
limited",  then  Ey  ■  0  at  y  ■  0  and  Cx  ■  0.  Then,  the  field  and 
Ion  density  profiles  are 
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and 


(19) 


(20) 


The  above  equation  for  oc  yields  an  Infinite  Ion  density 
at  the  mldplane,  y  ■  0,  as  a  result  of  the  neglect  of  the  Ion 
diffusion  term.  Ionic  diffusion  would  tend  to  decrease  the 
Ion  density  near  y  ■  0. 


According  to  Bq.  (19)  the  electric  field  strength  Is  zero 
at  the  plane  of  ionization,  y  ■  0,  and  increases  with  y1'  . 


When  the  Ionization  plane  is  a  corona  discharge,  another 
possible  boundary  condition  Is  that  Ey  ■  Ec  at  y  ■  0,  where  Ec 
Is  a  constant  equal  to  the  breakdown  rield  strength  of  the  gas 
(Ref.  5,  p.260).  Then  ■  Ec2  and 


(21) 

(22) 


The  solution  to  the  velocity  field  will  be  determined 
first  using  the  field  and  ion  density  relation  of  Eqs.  (19) 
and  (20).  From  Eq.  (13).  the  electric  field  component  in  the 
z  direction,  Ez,  is 


Substituting  Eqs.  (23)  and  (20)  into  Eq.  (10)  results  in 
the  differential  equation  for  the  velocity  w. 


uiiw  _  1  /S?  w  .  i£  -  0  .  (2K) 

9y  K  /  2K  y  82 


The  solution  of  Eq.  (2*1)  is  listed  in  Ref.  3  and  requires 
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a  numerical  evaluation  of  Integrals.  An  approximate  solution 
can  be  achieved  by  replacing  y1'2  by  Its  average  value  over 
OiySH,  which  Is  2H-1'2.  Equation  (2M)  then  becomes 


a2w 


1 

K 


_  IE  »  0 

SZ 


(25) 


This  Is  equivalent  to  taking  a  mean  value  of  the  charge 
density  distribution  represented  by  Eq.  (19)  which  yields 


and  then  substituting  Into  Eq.  (10). 

Using  H,  iowm,  wm  as  non-dimensionalizing  distance,  pres¬ 
sure,  and  velocity,  Eq.  (25)  In  non-dimensional  form  is 


afw  -  n  « -  !!se  ini .  0 

dn  Pc  2  d< 


(27) 


where  W,  p»,  n,  and  C  are  non-dimensional  quantities.  Np  is 
called  the  charge  number  and  is  equal  to 


°cm 
u 


mH  f2Jyc]  H2 

K  “L HK  J  UK  * 


v;mH 

The  Reynolds  number  NRE  is  defined  as  NREj  *  p 

u 


The  hydrodynamic  boundary  condition  is  W  ■  0  at  y  ■  *H, 

the  no  slip  condition  at  a  solid  wall,  and  ^[  «  0  at  y  ■  0, 

dn 

which  means  that  the  ionization  plane  imparts  no  shear  force  to 
the  moving  fluid. 


is 


The  solution  to  Eq.  (27)  with  the  above  boundary  conditions 


w  = 


-NRE  3P* 
2N0C  3  C 


1  - 


cosh  “pc _ n 


cosh  M 


pc 


(28) 


which  reduces  to  the  expected  parabolic  profile 
W  -  -  Nrf  _LE*(l-n2 ) 

3  C 


(29) 


\ 
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for  N  ■+  0. 

Pc 

Defining  a  friction  factor  f  ■  we  find  from  inte- 

3  C 

grating  Eq.  (28)  across  the  flow  area  that 


3/2 

f  .  2  2NpC _ _ 

NRE(NpC*  -  tanh  Npci)  » 

which  reduces  to  the  expected  f  ■  ^/Nre  for  plane  Poiseulle 
flow. 


(30) 


Velkoff  also  gets  an  approximate  solution  to  Eq.  (2*0 
by  taking  the  second  term  of  that  equation  equal  to  a  constant, 
its  average  value  over  0  1  y  1  H.  Equation  24  in  non-dimensional 
form  is  then 


32W  m  NRE  3p* 

jp-  -  "»c  -  —  rr  ■ 0  • 


(31) 


The  solution  of  Eq.  (31)  with  the  boundary  conditions  as 
before  is 


W 


npc  nRE 
2  4 


(32) 


The  friction  factor  f  *  -  dp/dj  is 


f 


6 

Nre 


which  again  becomes 


f  =* 


_6_ 

rIRE 


(33) 

(34) 


for  the  limiting  case  of  Npc  =  0.  Equations  (30)  and  (33)  are  both  approximate 
solutions  of  the  friction  factor  derived  from  Eq.  (24).  Equation  (30)  utilized 
an  average  value  of  y"'/'  and  Eq.  (33)  utilized  an  average  value  of  wy"’A  for  the 
solution.  They  can  be  regarded  then  as  integral  solutions  of  the  friction 
factor.  Figure  1  compares  the  graphs  of  the  two  integral  solutions  and  it  is 
seen  that  they  differ  by  less  than  11*,  up  to  a  charge  number  of  4. 


The  solution  to  the  velocity  field  will  next  be  determined  using  the  field 
and  ion  density  profile  according  to  Eqs.  (21)  and  (22).  Substituting  Eqs.  (22) 
and  (23)  into  Eq.  (10)  results  in  the  differential  equation  for  the  velocity  w. 
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dSf  vJy 

dy2  '  K2 


K« 


2Jyy  +  EC2K« 


(35) 


The  relative  magnitudes  of  the  two  terms  in  the  denominator  of  the  second 
term  will  be  determined.  For 

Jy  *  3*33  x  10"°  amp  in."2 

y  **  0.312  in.  max 

K  ■  0.216  in.2  volt"1sec"1 

€  »  2.5k  x  10-14  coul2lb-1in."2 

Ec  »  76.2  x  103  volt  in.-1 


Then 


and 


2Jyy  »  2.08  x  10"6  amp  in 


-1 


E2Kc  *=  2.84  x  10"4  amp  in.-1 
c 

Therefore,  since  2Jyy  is  negligible  compared  to  E2Ke,  it  will  be  neglected. 
Equation  (35)  then  reduces  to 
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wj 


“  Fe^  '  sf  =  0  ‘ 


(3o) 


In  nondimensional  form 


dfW 

3tj2 


JyH2 


2  bi 


(37) 


The  quantity  JyH2/nK2Ec  -  Ni  appears  in  Eq,  (37)  in  the  same  position  as 
Npc  in  Eq.  (27).  Since 


Jyn  _  _1_  j  J^rH 
4K2Ec  ~  Ec  L  2Ke 


WPp 


the  solution  to  Eq,  (37)  is  the  came  as  Eq.  (28)  except  that  Npc  is  multiplied 
by 
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Fig.  1  -  Friction  factor  times  Reynolds*  number  versus  charge  number  for 
theoretical  models 


-he  solution  to  Eq.  (37)  is 

u  „ .  * 

2Ni  df 


> 


(38) 


[ 


1 


cosh  Ni'*4  T) 
cosh  Nia 


and  the  corresponding  friction  factor  is 


f 


2Ni 


3/2 


NRE(Ni'A-tanh  Ni'4) 


(39) 


One  significant  difference  between  the  velocity  profile  solution  using  the 
ion  density  relation  of  Eq.  (26)  and  the  solution  using  (22)  is  that  in  the 
former  the  current  density  appears  in  the  non-dimensional  charge  number  to  the 
half  power,  whereas  in  the  latter  the  current  density  appears  to  the  fiyst 
power. 


The  assumption  is  sometimes  made  when  dealing  with  distribution  of  charges 
in  an  electric  field  that  the  charges  do  not  affect  the  electric  field.  Making 
this  assumption  then,  the  electric  field  strength  is  constant, 

Ey  =  H  0*0) 


where  V  is  the  voltage  applied  to  the  ionization  plane,  and  the  walls  are  at 
zero  potential. 

Now  from  Eq.  (12), 


Jy  ~  PcKEy 

P  =  JJL-  a  M  (41) 

KEy  KV 


Substituting  Eqs.  (40)  and  (23)  into  Eq,  (10)  results  in  the  differential 
equation  for  the  velocity  field  based  on  no  field  distortion  by  ions. 


d 


ap-., 

v 


KaV 


w 


(42) 


In  non-dimensional  terms,  Eq.  (42)  becomes 


d^W 

aq2 


JyH3 

nre 

fP 

p  w  - 
hk2v 

2 

26 

(43) 
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The  non-dimensional  parameter  Jy.H3/nK2V  »  N2  appears  in  Eq.  (43)  exactl; 
as  Npc  appears  in  Eq.  (27).  The  solution  therefore,  is  identical  to  Eq.  (28 
with  Npc  replaced  by  N2.  Since 


jyH3  1  r  jvh3  Y 

Mirv  "  V  L  2Ke  j  Lp°  * 


The  non-dimensional  parameter  in  the  solution  of  Eq.  (43)  is  Npc  multiplied  by 


if  act1  . 

V  L  2Ke  J 


The  corresponding  friction  factor  is 


_ _ 

Nre(N2,a  -  tanh  N2'A) 


As  in  Eq.  (37),  the  current  density  appears  in  the  non-dimensional  charge 
number  to  the  first  power. 

The  approximate  theories  discussed  are  tabulated  below: 


Charge  Density 
Distrubution 


[mT4 

L  KH  J 


Friction  Factor 


NRE(Npc''4  -  tanh  Npc”) 


Nop  V  3  ' 


Space  charge  limited  field 
at  midplane,  integral  solu¬ 
tion  of  Eq.  (24)  by  taking 
average  value  of  ym'A 


Space  charge  limited  field 
at  midplane,  integral  solu¬ 
tion  of  Eq,  (24)  by  taking 
average  of  wy",/l 


NnP(Niv*  -  tanh  Ni*) 


Critical  field  strength 
at  midplane 


Nre(N2I/1  -  tanh  N2*) 


No  distortion  of  field  by 
charges 
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where 


In  the  table  above  there  are  only  two  different  functional  relations 
between  friction  factor  and  charge  number,  which  are  shown  plotted  in  Fig.  1. 

If  f  =  dp/d?  from  Eq.  (30)  is  substituted  in  Eq.  (28),  there  results  the 
expression  for  the  non-dimensional  velocity  distribution  as  a  function  of 
charge  number. 


cosh  Npc'/*i} 
cosh  NPc4 


1 


tanh  Npc'/i 
NPC'A 


(1*5) 


Figure  2  shows  the  shape  of  the  velocity  distribution  from  Eq.  (45)  with 
charge  number  as  a  parameter.  The  velocity  distribution  becomes  flattened  in 
the  central  region  of  the  flow  channel  with  increasing  charge  number,  in  exact 
analogy  with  Hartmann  flow  in  magnetohydrodynamics. 


Kinetic  Theoi 


sroach 


The  kinetic  theory  approach  treats  the  ions  as  being  discrete  particles 
dispersed  within  a  fluid.  If  the  ions  have  a  statistically  different  motion 
than  the  molecules  of  the  main  fluid,  then  forces  on  the  main  fluid  arise 
because  of  collision  between  the  ions  and  neutral  particle. 


Considering  flow  between  infinite  parallel  plates  and  using  the  same 
notation  as  before,  the  equation  for  the  velocity  is,  from  Refs.  6  and  7 


H  rrr  =  X  +  -  f 2m*n kT  ] '/‘  dio2n0pc(w-wi)  (46) 

ay  oz  3 


where 
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m*  ■  reduced  mass  m^mQ/m^+mQ 
mi,  mo  ■  mass  of  ion  and  mass  of  neutral  particle 
T  «  absolute  temperature 
dio  -  1/2  (di+do) 

di»  dQ  *  diameter  of  ion  and  diameter  of  neutral  particle 
Pc,  n0  ■  number  density  of  ions  and  neutral  particles 
wi  =  statistical  average  velocity  of  ion 


The  last  term  in  Eq.  (46)  is  analogous  to  a  body  force  per  unit  volume  in 
the  hydrodynamic  formulation.  The  solution  of  Eq.  (46)  requires  a  specification 
of  wi,  the  ion  velocity.  The  greatest  contribution  of  the  last  term  of  Eq. 

(46)  occurs  for  Wi  =  0,  which  would  mean  that  the  ions  travel  across  the  flow 
with  no  downstream  motion.  This  is  exactly  the  same  simplifying  assumption 
used  earlier  when  Jz  was  taken  equal  to  zero.  The  solution  to  Eq.  (46)  will 
be  obtained  using  Wi  =  0;  the  friction  factor  thus  obtained  will  be  an  upper 
limit  on  the  effect  of  ionic  coupling  to  the  mean  f]ow. 

First  nondimens ionali zing  Eq.  (46), 

&  -  -  [2jTm*kTP  aio2  n0pci!n-!ll  i  =  0  .  (47) 

o'T  3  M  2 


Comparing  Eq.  (47)  with  Eq.  (27),  the  quantity 

|  [2rtm*kTP  dio2  n0  pc  -  K3  (48) 


is  exactly  analogous  to  Np  . 

Since  pc>  the  ion  number  density  can  possibly  be  a  function  of  y,  the 
solution  of  Eq.  (47)  requires  a  separate  solution  of  the  electrodynamic  equa¬ 
tions  to  obtain  pc.  This  requirement  is  the  same  as  the  requirement  of  obtain¬ 
ing  pc  in  the  hydrodynamic  formulation. 

For  the  assumption  of  pc  =  constant  =  Jy/KEC,  (which  led  to  nondimens ional 
parameter  Ni  in  the  hydrodynamic  formulation)  and 

mp  =  my  =1.02  x  10"25  LBM/Molecule 

m*  =  5.13  x  10-2C  LBM/Molecule 
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d^o  ■  di  ■  do  ■  1.48  x  10-®  in. 

T  -  300#K 

nD  ■  4.03  x  IO20  Molecules  in. 

U  ■  1.23  x  10-5  lbm  ft_1sec_1 
H  -  0.312  in. 

K  ■  0.2l£  in.2  volt“1sec“1 
Ec  **  76.2  x  103  volt  in.-1 

we  find 

N3  »  0.20  Ni  (49) 


and  for 


hi 

KH  J 


'A 


which  led  to  the  nondimen  5 ion&l  parameter  Np_  in  the  hydrodynamic  formulation, 
we  find 


N3  =  0.14  Npc  .  (50) 

The  kinetic  approach  leads  to  the  identical  problem  as  the  hydrodynamic 
formulation,  i.e.,  the  ion  number  density  pc  must  be  specified  for  a  deter¬ 
mination  of  the  parameter  equivalent  to  the  charge  number.  For  numerical 
values  of  the  quantities  resulting  from  the  kinetic  approach,  the  resulting 
nondimens ional  parameter  is  seen  to  be  less  than  the  analogous  nondimens ional 
parameter  resulting  from  the  hydrodynamic  approach,  based  on  the  same  ion 
density. 

The  simple  kinetic  approach  assumes  that  momentum  is  transferred  by  col¬ 
lision  of  ions  and  neutral  particles,  treated  as  hard  spheres  which  exert 
forces  only  on  collision.  Since  changed  particles  can  exchange  momentum  with¬ 
out  actually  colliding  in  the  narrow  sense,  close  agreement  in  the  determina¬ 
tion  of  the  charge  number  between  the  simple  kinetic  approach,  and  the  hydro- 
dynamic  approach  is  not  expected. 
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SECTION  III 


EXPERIMENTAL 


Description  of  Flow  Channels 

The  production  of  ions  in  all  phases  of  the  test  work  was  achieved  by  the 
positive  corona  discharge  mechanism.  A  high  positive  voltage  was  applied  to 
thin  steel  wires,  with  a  grounded  electrode,  usually  the  duct  walls,  located 
nearby.  The  resulting  electric  field  is  much  more  intense  near  the  wires  than 
near  the  ground  electrode,  and  if  the  applied  voltage  is  above  a  critical  value, 
breakdown  of  the  air  in  the  high  field  region  near  the  wire  results  in  ioniza¬ 
tion  of  the  air. 

Within  the  region  near  the  wire  where  ionization  takes  place  these  exist 
both  negative  and  positive  ions;  the  negative  ions  or  electrons  travel  to  the 
wire  and  the  positive  ions  travel  out  of  the  ionization  region  and  drift  toward 
the  grounded  electrode.  Except  for  the  thin  region  near  the  wire  then,  the 
fluid  under  investigation  consists  of  air  containing  drifting  positive  ions. 

Since  the  applied  voltage  is  the  cause  of  the  ionization  and  is  also  the 
cause  of  the  positive  ion  drift,  it  is  clear  that  the  amount  of  ionization  and 
drift  velocity  are  not  independent  for  this  method  of  ion-flow  tests.  Varia¬ 
tion  in  drift  velocity  and  ion  current  could,  however,  be  achieved  by  using 
wires  of  different  diameters. 

The  first  rectangular  channel  to  be  constructed  under  the  test  program 
was  a  l/4  inch  5  t  inch  x  3  foot  long  duct  with  one  5 -inch  side  as  the  grounded 
electrode.  Wires  (0.004-inch  dia)  were  stretched  along  the  opposite  5-inch 
side  in  the  flow  direction.  The  wires  were  spaced  l/4-inch  apart  and  located 
tangent  to  the  wall  surface. 


Wirtt 

///////////  U 


TTTTTTJJTTT7 


Dielectric  Wall 


71ft  TTTTT7 
Grounded  Wall 
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The  flow  channel  with  no  voltage  on  the  wires  would  closely  approximate 
flow  through  a  rectangular  duct,  and  serve  as  a  definite  limiting  case  for 
the  study  of  the  ion  coupling  phenomenon. 

With  the  application  of  slowly  increasing  voltages  to  the  wires  in  the 
above  channel,  it  was  observed  that  a  constant  current  was  not  obtainable,  but 
rather  frequent  arcing  and  spark  discharges  occurred.  At  this  time  it  was 
postulated  that  the  proximity  of  the  dielectric  wall  (Plexiglas)  was  in  some 
way  affecting  the  comona  discharge  mechanism.  To  further  investigate  this 
occurrence,  some  test  work  was  performed  to  measure  the  effect  of  a  dielectric 
surface  near  a  corona  wire.  This  work  is  described  in  Appendix  III  where  it 
is  shown  that  the  proximity  of  a  surface  to  a  corona  discharge  greatly  affects 
the  corona  characteristics. 

Since  the  main  object  of  the  study  was  the  ion  flow  coupling  with  air 
flow,  another  rectangular  duct  5/8  inch  x  5  inches  x  10  feet  long  was  con¬ 
structed  as  shown  in  Figs.  3-7* 

The  corona  wires  were  stretched  lengthwise  along  the  midplane  of  the 
channel,  thereby  obviating  the  surface  effect  noticed  in  the  earlier  flow 
channel. 

Static  pressure  taps  were  located  along  the  centerline  of  the  upper  5-inch 
wall  at  one-foot  intervals  along  the  duct  length.  Traversing  total  pressure 
probes,  constructed  of  glass,  were  located  at  distances  of  5,  9,  and  9.8  feet 
from  the  duct  inlet  (hereafter  referred  to  as  Positions  A,  B,  and  C,  respec¬ 
tively).  Probes  at  Positions  A  and  C  were  located  in  the  center  of  the  duct 
in  the  transverse  direction  and  hence  traversed  the  duct  cross  section  midway 
between  two  corona  wires;  the  probe  at  Postion  B  was  directly  over  a  wire  and 
hence,  only  traversed  half  the  height  of  the  duct. 

The  output  from  the  high-voltage  supply  was  connected  to  ten  variable 
resistors  in  parallel  (Fig.  5)«  Each  resistor  was  in  turn  connected  to  one 
of  the  ten  corona  wires  in  the  duct;  then  with  a  given  voltage  setting  on  the 
high-voltage  supply,  the  voltage  drop  across  the  variable  resistors  enabled  a 
different  voltage  to  appear  on  each  cornona  wire  by  adjusting  the  resistor. 

The  voltage  adjustment  was  necessary  since  each  wire  has  a  slightly  different 
ampere-voltage  characteristic  (Fig  8),  and  it  was  desired  to  maintain  equal 
currents  in  each  wire.  For  all  tests  performed,  the  resistors  were  so  adjusted 
that  a  uniform  current  existed  in  each  wire.  From  Fig.  8  then  it  is  seen, 
from  the  variation  in  voltage  at  a  given  current,  that  the  voltages  appearing 
on  the  wires  could  differ  from  each  other  by  up  to  150  volts.  The  range  of 
voltage  applied  to  the  corona  wires  for  all  test  work  on  the  duct  was  5000  to 
5700  volts  (Fig.  8). 


Electrical  Characteristics  of  the  Test  Apparatus 

Upon  application  of  the  high  voltage  to  the  corona  wires,  a  uniform  blue 
corona  glow  was  observed  and  repeatable  ampere-volt  characteristics  were 
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Cross  sectional  drawing  of  the  flow  channel 
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Electrical  circuit  of  the  flow  channel  tast  apparatus 


Corona  Current  Per  Wire 


A  ■ 


Fig.  8  -  Corona  current  versus  volts  applied  to  wire  showing  the  variation 
in  characteristic  among  the  wires 
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achieved.  Figure  8  shows  the  range  of  characteristics  of  the  different  wires, 
and  Fig.  9  shows  the  ampere-volts  characteristic  for  ascending  the  descending 
voltages . 

Figure  10  shows  how  the  ampere-volt  characteristic  acts  with  a  flow  of 
air  in  the  duct.  Since  the  ampere-volt  curves  at  various  Reynolds  numbers  are 
approximately  identical,  it  can  be  concluded  that  the  mean  air  motion  has  no 
effect  on  the  corona  voltage -ampere  characteristic.  It  may  also  be  stated 
that  within  the  range  of  air  flow  rates  tested,  there  is  no  loss  of  ions  by 
being  blown  from  the  duct  exit,  as  this  phenomenon  would  cause  the  ampere-volt 
characteristic  of  the  electrical  circuit  to  be  dependent  on  the  Reynolds  number. 

The  plane  of  the  ten  wires  is  positioned  at  the  midplane  of  the  channel, 
and  consequently  an  equal  current  should  flow  to  each  grounded  wall  if  the  wire 
spacing  were  exact.  For  the  experimental  tests  performed,  it  was  possible  by 
positioning  of  the  wires  to  obtain  ion  currents  to  the  two  grounded  walls  which 
were  equal  to  within  5$. 

Since  the  approximate  theory  developed  earlier  led  to  the  condition  that 
the  transverse  current  density  Jy  *  constant,  Jy  could  be  obtained  by  dividing 
the  total  corona  current  by  the  total  exposed  area  of  Loth  grounded  walls. 

The  nondimens ional  charge  numbers  as  developed  earlier  by  means  of  various 
approximations  are  tabulated  below  for  this  flow  channel  at  several  values  of 
total  corona  current. 


Total  Corona  Current 
in  Test  Channel 

„  JyH2 

Ni  =  -*4r— 
mK2Ec 

.  Jyl 

N2  = 

pK‘ 

0 

0 

0 

0 

1  ma 

3.53 

0.076 

0.3^3 

2.5  ma 

5.58 

0.1837 

0.836 

4.0  ma 

7.06 

0.304 

1.31 

5.5  ma 

8.29 

0.415 

1.76 

It  was  found  that  the  corona  wires,  which  were  bright  and  lustrous  in 
appearance  when  new,  exhibited  a  light-brown  coating  on  their  surface  after 
an  extended  period  of  corona  discharge.  This  phenohenon,  postulated  to  be 
caused  by  the  corrosive  action  of  ozone,  caused  a  slight  change  in  the  volt- 
ampere  characteristics  of  the  wire.  As  shown  in  Fig.  11,  the  effect  is  to 
require  approximately  150  more  volts  on  the  old  wire  tc  obtain  the  same  current 
as  on  a  new  wire. 


24 


Corona  Current,  Microampere 
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Kilovolts  Applied  to  Wire 


Fig.  9  -  Corona  current  versus  volts  applied,  for  wire  No.  5, 
showing  the  different  characteristic  for  increasing 
and  decreasing  voltages 
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Corona  Current >  Microampere 


Kilovolts  Applied  to  V/ ire 


Fig.  10  -  Corona  current  versus  volts  applied,  for  wire  No 
showing  that  Reynolds  number  has  no  effect 
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was  determined  by  measuring  the  static  pressure  decrease  along  the  duct  in  the 
flow  direction.  It  can  be  seen  from  Figs.  12,  13,  14,  15,  and  l£,  that  the 
pressure  gradient  in  the  flow  direction  appears  as  a  constant,  which  supports 
the  assumption  that  the  flow  field  can  be  described  as  "fully  established,"  or 
one  dimensional.  Figure  17  shows  the  log  friction  factor  plotted  against  log 
Reynolds  number.  The  curve  for  zero  current  is  parallel  to  and  slightly  above 
the  line  f  =  6/Nr£  which  is  the  result  for  simple  laminar  flow  of  air  between 
infinite  parallel  plates.  The  amount  by  which  the  zero  current  line  is  above 
this  curve  can  be  attributed  to  the  side  wall  and  wire  influence  on  the  fric¬ 
tion,  and  experimental  error. 

For  Reynolds  numbers  below  about  650,  it  is  seen  on  Fig.  17  that  the 
curves  for  friction  factor  with  ion  flow  are  parallel  to  and  above  the  curve 
for  simple  air  flow.  This  indicates  that  the  ion  coupling  effect  behaves  as 
a  laminar  flow,  i.e.,  friction  factor  times  Reynolds  number  equals  a  constant 
for  a  certain  charge  number  Npc< 

Above  Reynolds  numbers  of  about  650,  the  curve  in  Fig.  17  for  zero  current 
departs  from  linearity,  indicating  transition  to  turbulence.  The  flow  regime 
above  =  650  will  therefore  be  considered  the  turbulent  regime.  It  is  noted 
from  Fig.  17  that  within  the  turbulent  regime  the  measured  friction  factors 
with  ion  current  fall  on  the  same  curve  as  the  flow  with  no  ion  current.  In 
this  flow  regime  the  frictional  resistance  caused  by  turbulent  mixing  becomes 
large  in  relation  to  the  induced  electric  body  force,  and  thus  would  tend  to 
overshadow  the  measureable  effect  of  the  charge  number.  In  addition,  the  eddy 
diffusivity  associated  with  turbulence  might  sufficiently  affect  the  ion  dis¬ 
tribution  across  the  duct  to  result  in  smaller  electric  body  force  than  would 
exist  in  laminar  flow  at  the  same  charge  number. 

The  quantitative  effect  of  ion  current  on  the  friction  factor  is  depicted 
in  Figs.  18  and  19,  which  include  only  data  in  the  laminar  flow  regime.  Also 
included  on  Fig.  18  are  curves  representing  the  approximate  solutions  developed 
earlier  (Eqs.  (30,  33j  39>  44).  A.1 1  the  theoretical  curves  have  been  shifted 
upward  on  Figs.  18  and  19  by  a  facotr  of  1.53  20  that  they  agree  with  the 
experimental  friction  factor  data  at  zero  ion  current.  This  upward  shift 
represents  the  friction  contributed  by  the  side  walls  01'  the  channel  and  the 
corona  wires —  not  included  on  the  theoretical  treatment,  as  well  as  systematic 
experimental  errors. 

Figure  18  shows  that  the  experimental  data  lie  between  the  first  two  and 
last  two  of  the  approximate  solutions  listed  on  page  11.  The  upper  two  curves, 
Eqs.  (30)  and  (33)  are  functions  of  Npc,  which  contain  the  ion  current  to  the 
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Nondimensional  static  pressure  versus  distance  along  duct  for  a 
Reynolds  number  of  approximately  lUO,  with  charge  number  as  parameter 
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Fig.  15  -  Nondimens ional  static  pressure  versus  distance  along  duct  for  a 

Reynolds  number  of  approximately  540,  with  charge  number  as  parameter 
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Fig.  17  -  Log  friction  factor  versus  log  Reynolds  number, 
with  charge  number  as  parameter 
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Friction  factor  times  Reynolds  number  versus  charge  number,  showing  comparison 
between  experiment  and  a  theoretical,  model 


mt* 


one-half  power,  and  hence  appear  as  curved  lines.  The  lower  curves,  Eqs.  (39) 
and  (44)  are  functions  of  Ni  and  N2  which  contain  the  ion  current  to  the  first 
power,  and  hence  appear  as  approximately  straight  lines.  Although  Eq.  (44)  is 
closest  to  the  data  points  at  high  values  of  current,  it  is  to  be  noted  that 
the  data  functionally  does  not  behave  as  a  straight  line  for  small  values  of 
current . 

On  Fig.  19  the  experimental  friction  factor  data  is  plotted  against  Npc. 
It  is  seen  that  the  data  fall  fairly  well  on  a  straight  line  for  all  values  of 
charge  number,  and  hence  it  is  concluded  that  the  friction  factor  is  function¬ 
ally  dependent  on  the  ion  current  to  the  one-half  power.  Although  the  theo¬ 
retical  curves  on  Fig.  19  indicate  much  too  great  on  ion  effect  than  that 
measured,  it  must  be  remembered  that  they  are  dependent  on  the  value  of  the 
mobility,  which  appears  in  the  charge  number  to  the  3/2  power,  and  the  exact 
value  of  mobility  to  be  used  in  any  one  experiment  is  somewhat  uncertain.5  A 
straight  line  was  drawn  through  the  data  points  on  Fig.  19  which  represents 


fNRE  =  9.2  (1  +  0.13  Npc),  (51) 


an  empirical  line  which  best  fits  the  experimental  data. 


Effect  of  Ions  on  Transverse  Pressure  Gradient 


For  fully  established  flow  with  no  ion  effect,  the  static  pressure  trans¬ 
verse  to  the  flow  is  constant.  However,  in  the  presence  of  a  field  E  and  a 
charge  density  pc,  there  results  a  body  force  in  the  transverse  (y)  direction 
equal  to  pcEy,  which,  for  one  dimensional  flow  is  equal  to  the  transverse 
pressure  gradient,  bp/ by. 

Since 


dj  bjz 

from  Eq.  (6)  and  J_  was  assumed  zero,  there  results 


Jy  =  PcKEy  =  constant 


(52) 


and  since  pcEy  =  dp/dy  from  Eq.  (ll)  then 

a  'll.  =  constant. 
by  K 


(53) 
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Static  pressure  measurements  were  made  using  the  traversing  pressure 
probes  with  zero  air  flow  in  the  duct.  Figure  20  shows  the  pressure  varia¬ 
tion  at  locations  A  and  C  which  represent  traverses  between  corona  wires.  It 
is  seen  from  the  slopes  of  the  curves  that  the  experimen tally  measured  pres¬ 
sure  gradient  is  approximately  equal  to  the  theoretical,  gradient  except  for  a 
region  near  the  wires  representing  about  one-fourth  of  the  duct  volume.  In 
this  region  the  fact  that  there  are  discrete  wires,  each  acting  as  an  ion 
emitting  plume,  may  be  affecting  the  pressure  variation.  Therefore  it  can  be 
expected  that  the  transverse  pressure  gradient  near  the  plane  of  the  wires  will 
be  somewhat  different  from  the  gradient  far  from  the  wires.  It  is  also  noted 
that  the  pressure  distribution  at  location  C,  near  the  duct  exit,  is  similar 
to  the  distribution  at  location  A,  as  expected  since  the  probe  locations  only 
differ  in  tie  value  of  axial  position  z. 

It  was  thought  that  the  transverse  pressure  gradients  existing  in  the 
channel  could  result  in  a  circulating  flow  around  the  corona  wires,  (corona 
wind)  which  would  also  have  the  effect  of  increasing  the  friction  factor  if 
the  circulation  were  active  with  an  axial  flow.  Therefore,  transverse  pres¬ 
sures  with  no  net  axial  flow  were  measured  with  the  pressure  probe  oriented 
along  the  axis  (the  total  pressure  probes  used  in  the  flow  tests)  and  also 
with  a  probe  oriented  perpendicular  to  the  bottom  duct  wall.  If  there  was 
any  air  movement  in  the  duct,  there  would  be  a  different  pressure  distribution 
read  by  the  probes  in  the  two  orientations.  Figure  21  shows  a  comparison  of 
the  pressure  measurements  made  at  location  A,  which  shows  some  scatter,  but 
no  definite  differences  in  the  pressure  distribution  as  measured  by  probes  at 
the  two  orientations.  For  all  the  pressure  measurements  at  locations  A  and 
C,  the  minimum  pressure  was  at  q  =  ±  0.2  to  ±  0.3,  this  fact  not  explicable 
by  a  model  of  a  planar  ion  source  may  be  caused  by  the  effect  of  the  individual 
wires . 


Figure  22  shows  the  pressure  distributions  measured  at  location  B,  which 
is  located  directly  above  a  wire.  Again  there  is  some  scatter  in  the  measure¬ 
ments,  and  it  is  not  evident  that  the  probe  orientation  is  important  here. 

The  data  points  indicate  that  the  pressure  distribution  does  not  follow  the 
slope  of  the  theoretical  dp/3y  =  j/K  line.  It  is  suspected  that  here,  when 
the  probe  is  directly  over  the  wire,  the  presence  of  the  probe  affects  the 
corona  discharge  and  hence  affects  the  ion  density  i’'  the  vicinity  of  the 
probe.  Figure  23  shows  the  effect  of  current  on  the  transverse  electric  pres¬ 
sure  difference  between  the  point  q  =  -  0.5,  which  is  halfway  between  the  mid¬ 
plane  and  the  lower  wall,  and  the  pressure  at  the  upper  wall,  wher^  q  =  1. 

The  measurements  show  an  approximately  linear  increase  of  pressure  with  current 
density,  as  expected,  and  the  slope  of  the  data  points  is  approximately  equal 
to  the  theoretical  line  based  on  dp/py  =  Jy/K. 


Effect  of  lone  on  Velocity  Distribution 

Local  velocities  across  the  air  duct  were  determined  at  locations  A,  B, 
and  C,  by  measuring  the  difference  between  the  toted,  pressure  and  the  static 
pressure.  The  total  pressure  was  measured  with  a  glass  probe  oriented  up¬ 
stream.  The  static  pressure  at  the  upper  wall  at  locations  A,  B,  and  C  were 
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Fig.  23  -  Difference  'between  pressure  at  wall  and  probe  pressure  at  tj 
versus  corona  current  density 


obtained  from  the  pressure  at  a  point  away  from  the  wall  was  then  determined 
from  the  transverse  electric  pressure  curves  of  Figs.  20  and  22,  Since  no 
correction  was  made  for  the  effects  of  sheared  flow  on  the  static  probe  reading 
near  the  wall  or  wire,  errors  can  be  expected  for  the  velocity  determination 
close  to  the  wall  or  wire;  also,  because  of  the  large  scatter  in  the  trans¬ 
verse  pressure  gradient  measurements,  the  velocity  profile  data  are  subject  to 
a  large  experimental  error. 

Figure  24  shows  the  results  of  velocity  determination  for  location  A 
compared  with  the  velocity  profile  determined  by  the  approximate  theory.  The 
mean  velocity  wB  used  to  nondimens ionalize  the  velocity  was  determined  by  a 
measurement  of  the  total  flow  in  the  duct.  The  integrated  value  of  the  veloc¬ 
ity  curve  over  the  flow  cross  section  should,  therefore,  equal  1.0  if  no  exper¬ 
imental  erros  existed.  The  effect  noticed  is  that  the  experimentally  measured 
velocity  profile  is  flattened  much  more  than  the  theoretical  model  would  indi¬ 
cate. 


The  velocity  curves  for  low  Reynolds  numbers,  Figs.  25,  26,  27  in  the 
laminar  regime,  show  that  the  ion  coupling  acts  to  depress  the  velocity  at  the 
center  of  the  duct  and  causes  a  higher  velocity  gradient  (and  more  friction) 
near  the  walls,  in  qualitative  agreement  with  the  approximate  theory.  Near 
the  center  of  the  duct  the  deviation  of  the  theoretical  velocity  curve  and  the 
measured  values  is  the  greatest,  but  it  is  here  that  the  effect  of  the  discrete 
corona  wires  affect  the  model,  as  discussed  earlier  in  the  treatment  of  the 
transverse  pressure  gradient. 

The  velocity  curves  for  a  flow  rate  in  the  turbulent  regime,  Fig.  27 
shows  that  the  ion  effect  is  very  small,  in  agreement  with  the  effect  noticed 
on  the  friction  factor,  that  the  ion  effect  was  small  in  the  turbulent  regime. 

Figure  28  is  the  measured  velocity  distribution  at  location  B,  which 
shows  an  almost  flat  profile.  As  discussed  earlier  in  connection  with  the 
transverse  electric  pressure  at  location  B,  the  electric  pressure  determination 
here  is  in  doubt  because  of  the  probe  effect  on  the  local  corona  discharge. 

The  velocity  determination  at  location  B  is  therefore  thought  to  be  much  more 
in  error  than  at  locations  A  and  C  where  the  probes  are  between  wires.  Figure 
29  is  a  measured  velocity  distribution  at  location  C,  near  the  end  of  the  duct, 
and  is  seen  to  be  similar  to  the  results  obtained  at  location  A, 


Effect  of  Flow  on  Current  Density 

From  the  electrical  field  equations, 


V  x  E  =  0 


(54) 


and  in  this  channel  case  this  equation  reduces  to  Eq.  (7) 
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Fig.  25  -  Nondimens ionad.  velocity  profiles  for  Reynolds  number  of  265,  for  various  charge  numbers 
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(23) 

(55) 

(56) 

(57) 

(58) 

(59) 


This  above  analysis  indicates  that  unless  w  is  a  constant  value  across 
the  channel,  we  can  expect  that  a  variation  in  Ey  will  occur  along  the  axis 
of  the  channel.  If  one  examines  the  characteristics  of  a  corona  discharge, 
see  for  example  Fig,  9>  it  is  readily  apparent  that  the  discharge  is  highly 
sensitive  to  changes  in  voltage.  In  the  cirtical  discharge  region,  a  small 
change  in  voltage  will  lead  to  a  large  change  in  current.  Since  the  internal 
field  in  the  transverse  or  y  direction  is  related  closely  to  the  voltage,  it 
can  be  expected  that  if  changes  in  the  transverse  field,  Ey,  occur  along  the 
channel,  then  changes  in  the  local  transverse  current,  Jy,  should  result. 


50 


From  the  foregoing  analysis,  one  should  expect  to  see  an  approximately 
linear  change  in  Ey  along  the  channel.  If  this  occurs,  then  one  would  expect 
to  find  that  the  corona  currents  increased  along  the  length  of  the  channel. 
Uhder  such  conditions  for  a  very  long  channel,  it  could  be  expected  that  the 
air  flowing  along  the  channel  would  induce  very  high  local  Ey  fields  at  the 
end  of  the  pipe,  perhaps  sufficient  to  lead  to  breakdown.  Since  a  variation 
of  EL.  along  the  pipe  was  predicted  from  the  analysis,  it  was  decided  to  con¬ 
duce  a  test  to  determine  the  variation  of  current  density  along  the  channel. 

In  the  initial  tests,  the  grounded  channel  walls  were  constructed  of 
single  pieces,  and  measurement  of  local  current  density  was  impossible.  The 
flow  channel  was  later  modified  by  covering  one  grounded  wall  with  an  insu¬ 
lating  covering  and  painting  it  with  transverse  strips  of  conductive  silver 
paint,  each  strip  being  separated  from  the  adjacent  one  by  a  thin  unpainted 
area  approximately  l/32-inch  wide.  The  current  to  each  strip  divided  by  the 
area  of  that  strip  then  gives  the  local  current  density  at  the  location  of 
that  strip.  (See  Fig.  30)  These  modifications  were  made  to  determine  the 
extent  to  which  ion  density  is  affected  by  the  air  flow. 

Figure  31  shows  the  current  density  profile  for  zero,  265,  and  860 
Reynolds  number.  The  curves  show  only  a  slight  shift  in  shape,  toward  a  higher 
density  downstream  with  air  flow,  as  expected.  The  curve  for  zero  air  flow 
is  not  flat  because  slight  changes  in  Geometry  along  the  duct  due  to  wire  sag 
and  duct  alinement  greatly  affect  local  current  density.  The  isolated  points 
on  Fig.  31,  showing  low  values  of  current  density,  are  caused  by  the  effect 
of  nylon  filaments  stretched  across  the  duct  to  support  the  corona  wires  and 
eliminate  wire  sag. 

Figure  32  shows  the  average  current  density  for  an  area  of  duct  wall  in 
the  upstream  half  of  the  duct  and  an  area  of  wall  in  the  downstream  half  of 
the  duct,  as  a  function  of  Reynolds  number.  It  is  seen  that  the  air  flow  acts 
to  force  a  higher  current  density  downstream  and  a  lower  current  density  up¬ 
stream,  The  amount  of  change  depends  on  the  Reynolds  number  below  about  U00, 
and  is  independent  of  the  Reynolds  number  above  UOO.  These  curves  indicate 
that  an  ion-fluid  coupling  does  exist  and  suggest  that  the  postulated  flow 
induced  body  force  may  possibly  exist. 


Flow  Visualization  Attempts 

Smoke  was  injected  into  the  channel  through  the  total  pressure  probes  to 
obtain  a  visualization  of  the  flow  field.  (The  side  walls  were  transparent 
Plexiglas).  Although  smoke  particles  are  inherently  charged,  their  mass  is 
large  compared  to  an  ion,  and  therefore,  they  would  drift  to  the  electrodes 
much  more  slowly  than  the  ion  drift,  and  hopefully  the  smoke  particle  motion 
would  be  affected  by  the  local  fluid  motion  in  the  duct. 

With  no  ion  current,  the  sm.»ke  stream  provided  a  clearly  defined  visu¬ 
alization  of  the  streamlines ,  but  upon  the  onset  of  ion  current,  the  smoke 
pattern  was  rapidly  diffused,  mixed  with  the  air,  and  precipitated,  resembling 
a  turbulent  motion  in  appearance.  Consequently,  no  characteristic  flow  pattern 
was  discemable  with  the  smoke  studies. 
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Fig.  30  -  Photograph  of  the  flow  channel  with  one  wall  removed,  showing 
the  segmented  electrode 
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Effect  of  Fringing  Field 

Although  the  static  pressure  drop  along  the  flow  direction  was  meas¬ 
ured  as  a  straight  line,  indicating  that  the  ion  effect  was  continuously  active 
throughout  the  duct  length  and  was  not  an  end  effect,  the  following  test  was 
performed  to  determine  the  influence  of  the  electrical  end  effect,  or  "f ringing 
field."  Several  2-inch-wide  transverse  conducting  strips  at  different  locations 
along  the  length  were  disconnected  from  ground  and  left  "floating"  with  no 
current.  Air  passing  completely  through  the  duct  then  would  have  the  effect 
of  passing  through  six  channels  in  series,  as  far  as  the  electric  field  on  the 
lower  half  of  the  duct  is  concerned.  The  upper  grounded  wall  was  not  provided 
with  separate  conducting  strips  and  hence  was  at  ground  potential  completely. 

With  the  above  arrangement,  the  discontinuous  electric  field  had 
practically  no  effect  on  the  friction  factor,  (Fig  33)  and  hence  it  was  con¬ 
cluded  that  the  ion  coupling  effect  is  definitely  not  a  surface  or  "end"  effect. 

Effect  of  Number  of  Wires 


Throughout  all  of  the  test  work,  there  were  10  corona  wires,  spaced 
at  l/2-inch  intervals  across  the  5 -inch  width  of  the  duct*  To  determine  the 
importance  of  the  number  of  ion  producing  sites,  one  experimental  measurement 
of  friction  factor  was  done  with  only  five  wires,  spaced  at  one-inch  intervals, 
impressed  with  voltage.  The  other  five  wires  were  left  in  place,  but  discon¬ 
nected.  Figure  3^  shows  the  pressure  drop  along  the  duct  compared  with  a  test 
using  ten  wires,  both  tests  being  at  the  same  total  current.  The  corresponding 
friction  factor  for  the  five  wire  run  was  2l£  larger  than  the  ten-wire  case. 
The  difference  is  attributed  to  the  fact  that  with  one-inch  wire  spacing  the 
ratio  of  wire  spacing  to  half  height  of  the  duct  becomes  too  large  to  maintain 
geometrical!  similarity  with  the  ten-wire  case. 


Nondimensional  Axial  Distance,  Z/L 

•  Static  pressure  on  upper  vail  versus  axial 
distance  Z/L  shoving  effect  of  discontinuous 
electric  field  along  the  duct 
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Fig.  34  -  Static  pressure  on  upper  wall  versus  axial 
distance  Z  showing  the  effect  of  number  of 
corona  wires 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  effects  of  transverse  ion  flow  on  the  flow  of  air  in  &  rectangular 
duct  were  shown  to  display  the  same  trends  as  found  earlier3  in  experiments  on 
a  circular  pipe.  Friction  factors  were  significantly  increased,  velocity 
distributions  distorted,  and  a  useful  correlating  parameter,  the  charge  number, 
could  be  defined. 

Experimentally ,  it  was  found  that  the  laminar  flow  friction  factor  was  an 
almost  linear  function  of  the  nondimens ional  charge  number,  Npc,  as  predicted 
by  one  of  the  approximate  solutions.  The  data  was  well  represented  by  the 
relation 


f  «  2jI  (i  +  0.13  Npc) 

nRE 


up  to  a  charge  number  of  8.2,  which  was  the  greatest  charge  number  investigated. 

The  transverse  ion  flow  was  found  to  generate  an  electric  pressure  gra¬ 
dient  transverse  to  the  air  flow  which  agreed  closely  with  the  theoretical 
expression  for  a  planar  ion  source, 

3y  K  * 

near  the  channel  walls.  Near  the  midplane  of  the  channel  where  the  corona 
wires  were  located,  the  above  expression  was  not  applicable  because  of  the 
individual  effect  of  the  wires  which  deviated  from  the  ideal  planar  ion  source. 

Velocity  distributions  across  the  flow  channel  were  significantly  dis¬ 
torted  by  the  ion  current.  Typically,  the  centerline  velocity  was  decreased 
and  the  velocity  gradient  at  the  walls  increased. 

Local  transverse  ion  ’••urrent  density  at  the  channel  wall  was  measured 
and  found  to  be  dependent  on  the  air  velocity  within  the  channel.  The  effect 
of  the  air  flow  was  to  cause  an  increase  in  current  density  in  downstream 
parts  of  the  channel  and  a  decrease  in  upstream  parts.  At  sufficiently  high 
velocity  no  further  change  in  current  density  was  noted. 

It  is  recommended  that  further  work  be  done  to  determine  the  type  of 
possible  ion-fluid  coupling  mechanisms  and  the  extent  to  which  they  are  active 
in  specific  fluid-flow  configurations.  Included  among  the  various  possible 
ion-fluid  coupling  mechanisms  we  find  (l)  the  direct  generation  of  electric 
wind,  causing  mass  motion  along  the  direction  of  the  electric  field,8  (2)  the 
generation  of  secondary  flows  induced  by  electric  body  forces,  (3)  the  modi¬ 
fication  of  the  transport  properties  of  the  fluid  medium  caused  by  *.on  motion, 
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(4)  the  generation  of  flow- Induced  electric  body  forces,  which  Is  treated  in 
this  yeport  and  elsewhere,2*9*8  the  destabilization  of  a  fluid  interface  when 
two  fluids  or  phases  are  present. 10 * 11  Such  mechanisms  are  distinct  from 
effects  caused  by  a  non-uniform  electric  field  in  fluids  with  a  dipole  moment.12 

Experimental  and  theoretical  work  is  in  progress  on  the  study  of  flow 
instabilities  caused  by  electric  body  forces,  item  (2)  above,  which  has  been 
discussed  in  the  literature  as  formation  of  toroidal  cells19  or  longitudinal 
roll  cells,19*14  and  its  analogy  with  magnetic  fields.19  Such  secondary  flows 
have  analogies  in  classical  fluid  mechanics  where  the  destabilizing  body  forces 
are  caused  by  centrifugal  force  fields  or  gravitation  force  fields*19 

It  is  anticipated  that  further  experimental  work  will  include  the  appli¬ 
cation  of  electrostatic  probes  to  the  direct  measurement  of  ion  density  and 
field  gradient  within  a  moving  fluid,  and  the  development  of  flow  visualization 
techniques  suitable  for  use  in  electric  fields. 
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APPENDIX  I 


HYDRODYNAMIC  AND  ELECTRODYNAMIC  EQUATIONS 


Momentum  Equation 

y  direction 


Mass  Continuity 


Charge  Continuity 


Body  Force 


Maxwell's  Equations 


z  direction  P  ~r  ■  “ 

Dt 


+  +  +r* 


P  ,  d(pu)  ,  d(pw)  „ 


dpC  bjy  bjz 


Fy  -  »c%  *  MJ*H)y  ♦  |  p  v[e.^  ^  )TJ, 

•ieX|)p| 

Fz  -  PCEZ  ♦  Pe(J*H)2  ♦  |  »  v\ E.^  ||  )JZ 


(66) 

(67) 
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Equation  of  State 


APPENDIX  II 


APPROXIMATE  ANALYSIS  SATISFYING  THE  IRROTATIONAL  ELECTRIC  FIELD  CONDITION 

By  formulating  the  electro-fluid  dynamic  equations  as  a  one-dimensional 
problem,  it  was  shown  previously  (Eq.  (9))  that  Eq.  (7),  which  is  the  condi¬ 
tion  that  curl  E  «  0,  leads  to  the  result  that  the  velocity  is  constant,  and 
therefore  zero  everywhere.  The  analysis  was  continued,  but  for  a  finite 
velocity  the  approximate  theory  violates  the  electrodynamic  Eq.  (7),  the  condi¬ 
tion  of  an  lrrotational  voltage  field. 

In  the  following  analysis,  Eq.  (7)  is  alway  satisfied,  i.e.,  curl  E  •  0. 
and  the  field  is  a  potential  field,  but  it  results  in  a  violation  of  Eq.  (iS), 
which  is  the  relation  governing  the  electric  field  distortion  by  space  charges. 

From  Eq.  (5) 


EZ  -  -  i  .  (23) 

K 

Substituting  in  Eq.  (7) 


dz  K  dy 


0 


(72) 


take 


Since  w  is  a  function  of  y  only,  Ey  must  be  linear  in  z.  Therefore,  we 


% 


f(y) 


(73) 


where  0  is  a  constant  and  f(y)  is  a  function  of  y. 

For  the  y  dependency  of  Ey,  we  take  the  relation  of  Eq.  (19)>  based  on 
space  charge  limited  current  and  zero  velocity.  Then 


(74) 


Substituting  into  Eq.  (72)  and  integrating, 
v.f  |  [  Hd  ]*  |>  -  y3/z] 


(75) 
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“which  is  an  expression  for  the  velocity  distribution  obtained  without  consid¬ 
eration  of  the  fluid  momentum  (Eq.  (10) )  •  If  we  put  Eq.  (75)  into  Eq.  (23) 
and  Eq.  (74)  into  Eq.  (12)),  it  is  possible  to  obtain  an  expression  for  the 
electric  body  force  term,  pcEz. 


PcEz 


(76) 


Using  this  force  term  in  Eq.  (10)  and  integrating  Eq.  (10)  twice, 

*  -  fsfe  [yS '  8H3/2y3/2  +  ^  t  [y8  -  “*]  (77) 


where  the  first  term  is  the  electric-induced  term  and  the  second  term  is  the 
classical-parabolic  term. 

Calculating  the  mean  velocity  from  Eq.  (77) ,  the  friction  factor,  can  be 
calculated. 


f 


(78) 


where 

M  -  drift  number  -  [  1*  -  drlft  Teloclty  . 

>-  9EWn»2  J  mean  velocity 


If  Eq.  (74)  and  Eq,  (12)  are  combined  into  Eq.  (l6),  a  contradiction 
results;  therefore,  the  above  analysis  does  not  satisfy  the  condition  of  Eq. 
(16). 


The  above  analysis  indicates  that  the  friction  factor,  f,  varies  with  Jy 
and  that  fN^g  varies  with  It  was  verified  experimentally,  however,  that 

f  varies  with  J#  and  fNRE  is  constant  in  laminar  flow,  and,  therefore,  the 
above  analysis  is  not  a  realistic  model  of  the  ion-fluid  coupling. 


APPENDIX  III 


EFFECT  OF  NEARBY  SURFACES  ON  CORONA 

In  the  early  channel  flow  testa,  corona  wires  were  stretched  along  a 
plastic  surface  (Plexiglas)  which  formed  one  wall  of  the  duct.  It  was  found 
that  the  corona  current  characteristic  was  non-uniform  and  frequently  not 
repeatable.  Since  it  was  thought  that  small  variations  in  the  distance  between 
the  corona  wires  and  the  dielectric  surface  might  be  Important  in  the  corona 
current-voltage  characteristic,  an  apparatus  was  constructed  to  determine  the 
magnitude  of  the  effect  of  a  nearby  dielectric  surface  to  a  corona  discharge. 
The  apparatus,  shown  below,  maintained  a  fixed  distance  between  the  plane 
grounded  electrode  and  the  corona  wire  while  the  distance  from  the  corona  wire 
to  a  Plexiglas  surface  could  be  varied  by  a  screw  adjustment. 


A  0.010-inch  wide  x  l/32-inch-deep  slot  was  cut  in  the  Plexiglas  surface  to 
enable  the  surface  to  move  beyond  the  0.004-inch-diameter  wire.  Figure  35 
shows  the  effect  of  d,  the  distance  from  the  wire  to  the  plastic  surface.  It 
is  seen  that  for  small  values  of  d  the  current-voltage  curves  become  steeper, 
and  sparking  from  the  wires  to  the  ground  electrode  occurs  at  a  lower  voltage. 
For  d=0,  the  result  was  that  no  corona  current  occurred  at  all  voltages  up  to 
the  sparking  voltage.  Therefore,  the  effect  of  the  plastic  surface  was  very 
critical  to  the  corona  characteristic  for  small  values  of  d. 

It  was  desired  to  also  determine  the  effect  of  an  air  flow  over  a  corona 
discharge  while  being  affected  by  the  proximity  of  a  dielectric  surface.  For 
this  purpose  a  needle  utilized  as  a  corona  point  was  maintained  concentric 
within  a  hole  in  a  nylon  spool.  The  spool  could  be  raised  or  lowered  to  allow 
various  lengths  of  the  needle  to  protrude  beyond  the  hole  in  the  spool,  as 
shown  below.  The  distance  from  the  grounded  plane  electrode  to  the  needle  at 
positive  voltage  was  held  constant  while  the 
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distance  d  could  be  varied  and  the  air  flow  through  the  hole  in  the  spool  could 
be  metered.  Figure  36  shows  the  effect  of  distance  d  on  the  (ampere)  voltage 
characteristic,  which  shows  that  for  the  needle  tip  closer  to  the  nylon  spool 
the  current -voltage  curve  is  steeper  and  that  sparking  occurs  at  a  lower  volt¬ 
age  in  agreement  with  the  previous  results  on  the  wire. 

Figure  37  shows  the  effect  on  the  corona  characteristic  of  air  flow 
through  the  hole  in  the  spool.  As  indicated  on  the  figure,  at  the  end  of  the 
test  of  Fig.  37,  the  current  at  zero  air  flow  was  measured  again  and  found  to 
be  8.5u  amp  where  it  was  13.9h  amp  »t  the  beginning  of  the  test.  This  diffi¬ 
culty  of  a  changing  corona  characteristic  of  the  needle  point  is  a  result  of 
the  wearing  away  of  the  needle  tip  with  current,  which  changes  the  radius  of 
curvature  of  the  needle  tip.  A  change  in  curvature  then  changes  the  field 
strength  at  the  needle  tip  and  changes  the  current-voltage  characteristic. 

Because  of  the  effect  on  the  needle  tip,  repeatable  results  were  not 
obtainable.  However,  it  was  found  that  there  is  an  effect  of  the  air  flow  on 
the  corona  current.  Figure  37  shows  that  the  current  increased  with  higher 
velocity  air  flows,  porbably  resulting  from  a  flow-induced  change  in  the  space 
charge  around  the  needle.  It  is  to  be  noted,  however,  that  the  air  velocities 
of  Fig.  37  go  up  to  500  feet  per  second,  while  the  flow  velocities  used  in 
subsequent  channel  flow  tests  were  about  1  to  10  feet  per  second.  For  less 
then  10  feet  per  velocity,  the  slope  of  the  data  points  in  Fig.  37  would  indi¬ 
cate  a  negligible  charge  in  corona  current.  Indeed,  Fig.  6  shows  no  effect 
of  flow  on  the  corona  characteristic  for  Reynolds  numbers  of  up  to  1270, 
corresponding  to  a  flow  velocity  of  up  to  10  feet  per  second. 
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Air  Velocity  Through  Hole,  Feet  Per  Second 


Corona  current  versus  voltage  applied  to  a  needle 
point,  showing  the  effect  of  a  superimposed  air 
stream 
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Friction  factors,  velocity  distributions,  pressure  gradients,  and  density 
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channel.  Experimental  data  were  collected  to  determine  the  influence  of  the  ion 
flow  on  the  air  flow  and  the  extent  of  the  ion-fluid  coupling. 

The  electrodynamic  and  hydrodynamic  equations  were  set  up  and  several  approxi¬ 
mate  solutions  were  obtained  based  on  the  postulation  of  a  flow  induced  body  force 
and  the  assumption  of  fully  established,  one-dimensional  flow.  The  magnitude  of 
the  ion  effects  was  found  to  be  significantly  dependent  on  the  ion  density  distri¬ 
bution  within  the  flow  channel. 
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fied)  "The  Effect  of  Transverse  Ion  Current  on  the  Flow  of  Air  in  a  Flat  Duct," 
by  E.  Pejack  and  H.  Velkoff. 


Page  3 


Page  4 


Page  6 


Item  I  - 
Item  II  ■ 

Item  I  - 
Item  II  - 

Item  I  - 
Item  II  - 


Change  last  term  of  Eq.  (2)  to  p_Ey 

-  Filth  line  after  Eq.  (8)  change  dpOz  to  dp/dz 

Change  Ey  to  Ey 

-  Change  Jyy  to  Jyy 

Change  exponent  on  bracket  to  +  ^ 

-  Line  2  after  Eq.  (25)  change  Eq.  (19)  to  Eq.  (20) 


Item  III  -  Change  right  side'  of  Eq.  (26)  to  read 


/2Jy«W2 

\H 57 


Item  IV  -  Equation  following  Eq.  (27)  right  side  should  read 


f2Jyc]  A/2  j£ 

Lurj 


Item  V  - 


Page  7  Item  I  - 


Item  II 


In  Eq.  (28)  change  material 
inside  bracket  tc  read 


V2 

-  «<>•*»  Hpo  * 

co«h  N 


Line  following  Eq.  (32)  change  to  read  dp*/di 

N 

■  Equation  (33)  —  close  parentheses  after  •*- 


Item  III  -  Change  numerator  of  right  side  of  Eq.  (30)  to  read  2N_ 

Page  8  Item  I  -  Eq.  (37)  Change  tO  lpV.*t 

Page  10  Item  I  -  Eq.  (38)  Change  *P/P  C  to 

Item  II  -  Eq.  (43)  Change  »p/J  (  to  »p#/?  C  *»  «*i 


Page  11  Item  I  -  First  item  in  first  column  of  table  should  read 


Page  12  Item  I  -  Second  paragraph  change  f  ■  dp/d(  to  f  ■  -dp*/d{ 

Item  II  -  Eq.  (45)  put  exponent  “ *  on  second  set  of  brackets 


Page  13  Item  I  - 
Page  14  Item  I  - 
Page  15  Item  I  - 


Ordinate  should  read  tj  •  y/ti 
■  Change  t0 

Following  Eq.  (49)  change  Jy  to  “y 


D  D  C 

iPr^rmn  nra 

MAY  18 1967 
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ERRATA 

Page  16 
Page  28 
Page  30 
Page  38 
Page  45 


Page  46 
Page  66 


Page  70 


for  TECHNICAL  REPORT  #  2  December  16,  1966  Contract  DA-31-124-AR0-D-246 


Item  I  -  Fourth  Paragraph,  11. 


2,  delete  t  after  5 


Item  I  -  Delete  H  in  equation  at  top  of  page 


Item  I  -  Ordinate  should  read  p* 

Item  I  -  Third  paragraph,  last  line,  should  read  dp/dy 

ONfte  -  0 

Item  I  -  Label  symbols  . 

An,,  ■  T.08 
On,,  •  3.5J 

Item  II  -  Write  "Location  A"  on  graph 

Item  I  -  Change  Ngg  «=  530  to  read  =>  540 

Item  I  -  Change  Eq.  (77)  —  dp/dt  should  read  dp/dz 

Item  II  -  Change  Eq  (78)  —  N  should  read  N 

Item  III  -  In  equation  following  Eq,  (78)  change  E  to  €  and 


